Weak Coupling Phases of the Attractive t — t' Hubbard Model at the Van Hove Filling 
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Q\ , We apply a wilsonian renormalization group approach to the continuum limit of the attractive t—t' 

0^ ' Hubbard model, taken when the Fermi level is at the Van Hove singularity of the density of states. 

The model has well-defined scaling properties and the effective couplings display an unbounded flow 
J_h ■ in the infrared. We determine the leading instabilities by computing the different response functions 

up to i! = 0.5t. The phase diagram shows a large boundary between superconducting and charge- 
density-wave phases, that merge in a triple point with a phase separation instability. The latter is 
realized down to very low coupling constant, as the Fermi sea degenerates towards a pair of straight 

(N ■ 

lines near t' = Q.5t. 
71.27. +a, 74.20.Mn 

During recent years there has been important progress in understanding the properties of quantum electron liquids 
| | in dimension D < 3. One of the most fruitful approaches in this respect springs from the use of renormalization group 
(RG) methods, in which the different liquids are characterized by several fixed-points controlling the low-energy 
properties. The Landau theory of the Fermi liquid in dimension D > 1 can be taken as a paradigm of the success 
C$ ■ of this program. It has been shown that, at least in the continuum limit, a system with isotropic Fermi surface aijd 
regular interaction is susceptible of developping a fixed-point in which the interaction remains stable in the infraredul. 
This example is of particular significance since there are systems that do not lie in the perturbative regime, and yet 
we know that their description in terms of the effective Fermi liquid theory is correct, as its predictions apply fairly 

q ■ The question of whether different critical points may arise at dimension D = 2 is now a subject of debate@~@. The 
i— 1 1 interest has rosen in parallel to the failure of a global theoretical understanding of the high-T c superconductivity of 
the cuprates. The two-dimensional Hubbard model has been proposed long time ago as the correct starting point to 
k^J study the electronic correlations inside the copper oxide planeso, and several authors have also proposed recently that 
it should have a quantum critical point separating antiferromagnetic and d-wave superconducting phasesQ. However, 
' the elucidation of the existence of such critical point is quite difficult, specially in the strong coupling regime, since the 
presence of the strong correlations do not make obvious even the determination of the continuum limit of the model. 
C*~) ' Actually, apart from the above mentioned case of regular Fermi surface, in general it is not possible to find by simple 
inspection of the microscopic hamiltonian a set of scaling operators that behave properly under RG transformations. 
This is particularly true in the case of the Hubbard model at half-filling. 

The above considerations stress the interest of identifying microscopic models with correct scaling properties from 
the RG point of view. The t — t' Hubbard model filled up to the level of the Van Hove singularity provides such 
an example. The determination of the scaling properties and low-energy phases for repulsive interaction has been 
accomplished in Ref . || (see also Ref. |^) . We are interested now in the case of the t — t' Hubbard model at negative 
"-^J ■ U, which highlights some important properties of the electron interactions near a Van Hove singularity. 

We briefly review the way in which the low-energy_iimit is taken in the model, and how in this case a set of scaling 
O ■ operators can be read at once from the hamiltonianEB. High-energy and low-energy electron modes are separated by 
an energy cutoff E c , that is sent progressively towards the Fermi line as high-energy modes are integrated out in the 
RG process. When the Fermi level is at the Van Hove singularity, as shown in Fig. [j], most part of the low-energy 
states close to the Fermi line are concentrated around the saddle points at (it, 0) and (0, it), as these features are at 
the origin of the divergent density of states. Therefore, in building up the low-energy effective theory we may focus 
on two patches around the respective saddle points A and B, where the dispersion relation can be approximated by 

_ i 

£A,s(k) » T(t T 2t')kla 2 ±(t± 2t')k 2 y a 2 (1) 

a being the lattice constant. In the continuum limit the rest of modes are irrelevant, from the RG point of view. In 
fact, the effective action for the low-energy modes restricted to the region \e a (k)| < E c is given by 

S = J dud 2 k^2 ( w a ^ -( k ,^)<w(k,w) - e Q (k) a+ .(k,a;)a Q , iff (k,w)) 

0,£T 

-U / dud 2 k Pl (k,uj) Pl (-k,-uj) (2) 
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where a ottCr (a^ t a ) are electron annihilation (creation) operators {a labels the Van Hove point) and p^i are the density 
operators in momentum space. Under a change in the cutoff E c — > sE Cl with a corresponding scaling of the momenta 
k — > s x / 2 k, one can check that the effective action remains scale invariant after an appropriate scale transformation of 
the electron modes E3. We have therefore a microscopic model with a well-defined scaling behavior, 

susceptible of being studied by means of RG methods. 

As E c is sent towards the Fermi level and the electron modes are labelled according to the saddle point they 
are attached to, the interactions may be classified in the form shown in Fig. |[ The four possible types of local 
interaction L T intra, b T inter, Uback and U um k are renormalized by quantum corrections, following a conventional pattern 
from the quantum field theory point of view. By integrating high-energy excitations in the slices E c — dE c < \e\ < E c , 
the lowest order 0(dE c / E c ) corrections are given by the particle-hole diagram of Fig. ||(a). It is worthwhile to 
stress that, within this wilsonian RG approach, there are no more diagrams renormalizing the interaction coupling 
constants, since the particle-particle channel produces a contribution that is in general ~ (dE c ) 2 . This point has been 
conveniently clarified in Ref. j| It is only when the total momentum of the colliding particles equals zero that the 
particle-particle channel develops a contribution ~ dE c , but this points at the correction of a correlation function, 
defined at a particular momentum value, rather than at the renormalization of a coupling constant. 

As remarked in Ref. ||, the diagram in Fig. |^(a) corresponds to an antiscreening effect, that tends to enhance 
repulsive interactions. By the same reason, its effect in the present model is that of reducing the attractive interaction 
and cannot produce by itself any instability driving away from normal metallic behavior. The renormalized on-site 
interaction has the asymptotic behavior U ~ 1 / | log i? c | as the high-energy cutoff is reduced. As observed in Ref. 
|Ti| , however, there are effective interactions generated by quantum corrections that are not present in the original 
hamiltonian. These are interactions between currents with parallel spin, that appear through second order processes 
like that shown in Fig. |3|(b). They are strongly momentum dependent, as the polarizabilities near q = and 
q = Q = {tt, tt) are given respectively by 



x(q,w = 0) = 
x'(q,^ = o) 



2ttH 



2tt 2 < 



log 



log 



E 



ta 2 (q~Q) 2 



(3) 
(4) 



where c = l/^/l — 4{t' /t) 2 and c' = log (l + y/l — 4(i'/t) 2 J /(2t'/t) . These interactions give rise to a potential 

that is actually singular at small momentum transfer. In real space this corresponds to an interaction decaying like 
~ 1/r 2 . It is not strange that this kind of interaction may arise in the renormalization of the model since, together 
with the purely local interaction, it corresponds to the other not irrelevant four-fermion operator that may appear 
in the effective action. According to the wilsonian RG approach, once these effective interactions appear they have 
to be considered on the same footing than the original bare interactions, since they are needed for the complete 
renormalization of the model. 

The processes shown in Fig. |^(b) are singular at small momentum transfer and at momentum transfer ~ Q = {tt, tt). 
Therefore we have to introduce new couplings Vi n tra, V inter, Vback and V um k for the effective potential between electrons 
with parallel spins, according to the classification in Fig. |^. The important point is that these new interactions are 
also attractive since they are due to over screening, that is, to the screening of an interaction vanishing at the classical 
level. 

The overscreening effect can be also understood within the RG framework, by solving the flow equations with the 
initial conditions Vi n tra(E c ) = Vi n ter{E c ) = Vback{E c ) = V um .k{E c ) — at the upper value E c of the cutoff. The 
most significant contributions to the renormalization of the couplings are given by the diagram in Fig. |(b) at small 
momentum transfer, for Vintra and Vinter, and at momentum transfer Q, for Vback and V um k- In the diagram the 
interaction lines may stand either for the local interaction potential or for the 1/r 2 potential. One has to realize that, 
under renormalization, 1/r 2 effective interactions may also arise between currents with opposite spin, that we denote, 
conserving the previous notation, by V±i n tra, V±i n ter, V±b a ck and V± um k- The RG flow equations, that reflect the 
change of the couplings by integration of particle-hole processes at the high-energy cutoff E c , are given by 



Q 1 / 2 2 2\ 

n ~ {Vintra i Vinter) 7; T~ c ( {Uintra i Winter) ~t~ {Vintra i Vi n ter) ~t~ {V^intra i V\_i n ter) ) 

OE c 2lt l t V / 



(5) 



dE c 



2irH 



C (2 {Uintra i Uinter) {Vintra i Vinter*) ~t~ 2 {Vj_i n tra i Vj_i n ter) {Vintra i Vinter)) (6) 



and two more equations obtained from the former two by the replacements c <-> c', Xi n tra Xback and Xi n ter ^ X um k 
{X = U, V, Vj_)- Given that at the begining of the RG process all the interactions vanish, except Uintra, U inter, Uback 
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and U um k that equal the value of the local attraction U, it is easily seen that the instabilities of the model are given 
by the flow of the V and V± couplings, that become increasingly attractive at low energies. 

The instabilities of the coupling constants characterize the different ground states of the system, which may be 
determined by studying the different response functions. These are correlators computed at particular values of the 
momentum, which signal the way in which symmetry breakdown takes place in the model. One may check that in 
the present case of attractive interaction (U < 0), the only correlators that may diverge at a given frequency are 

those for the operators J2k ( a A|(k)aJ 1 |(-k) + a^ T (k)a^ (— k) + h.c\ pj(Q,a;) + pj(Q,w) and p|(0,cj) + p^(0,w). 

Their response functions characterize, respectively, s-wave superconductivity (i?scs), charge density wave (Rcdw) 
and phase separation (Rps) instabilities. 

We compute the response functions by exploiting again the scaling properties of the model. For this purpose we 
establish the dependence of the correlators on the cutoff E c and, taking into account the scale invariance of the 
model (up to logarithmic renormalizations of the couplings), we introduce the dimensionless scaling variable E c /lu to 
determine the dependence on the frequency to. The starting point is the perturbative computation of each response 
function, that exhibits a logarithmic dependence on E c /u>. In contrast to the RPA, we do not sum up the iteration of 
particle- hole bubbles, as thiSpGannot be a reliable expansion of the correlator in a model with strong renormalization 
of the one-particle propertiesfcj. Instead we compute the variation of the correlator under a reduction of E c and write 
it down in terms of the correlator and coupling constants at th&mew value of the cutoff. This procedure is similar to 
that followed in the study of one-dimensional electron systemal3 or coupled one-dimensional chainstil 

The scaling equations in the present case turn out to be 

^ W — TT~F TL (Vback + Vumk + V±back + V±_ um k) -f^RcDW (J) 

oE c n 2 t E c n 2 t E c 

dRscs _ c \og{E c /uj) c \og(E c /uj) 

~dE~T ~ ~2^t E~ c M2i<y±**tra + V ±umk )— ^ R S Cs (8) 

n ,-, — = n— (Vintra + Vinter + V±Antra + V^inter) ~^-RpS (9) 

oE c ix z t E c TT 2 t E c 

The equation for Rscs shows the nontrivial scaling factor log(E c /u>), as a consequence of the enhanced susceptibility 
at zero momentum in the particle-particle channel, that behaves as ~ log 2 (E c /u>). Such factor is nothing but a 
reflection of the divergent density of states at the Van Hove singularity, and it does not spoil, in any event, the scaling 
properties of the model. 

As remarked before, the natural scaling variable in the model is E c /u, what becomes now evident by inspection of 
the scaling equations (0)-(^|). The fact that the renormalized coupling constants follow an unbounded flow may cast 
doubts in the reliability of the perturbative RG approach. The divergences that, as a consequence of it, are found in 
the response functions at certain values of E c /lu provide, however, a way of discerning the competition between the 
different instabilities of the system. The response function having the strongest divergence as the value of u> / E c is 
lowered characterizes the ground state of the system. According to this criterion, we have determined the different 
phases at weak U coupling and t' up to 0.5t, where the Fermi sea degenerates into a pair of straight lines. The 
phase diagram is shown in Fig. ^ . We recall that the charge-densit^-wave and superconducting ground states are 
degenerated in the t' = attractive Hubbard model at half-fillingEio. On the other hand, it has been shown that 
phase separation cannot take place in the Hubbard model on a bipartite lattice, at any value of the interactionlij. 
The proof of this statement does not follow in our case, however, as the next-to-nearest neighbor hopping t' spoils 
the splitting into two sublattices. 

Our results should be applicable in the weak coupling regime and to Fermi surfaces corresponding to an appreciable 
value of t' . This is because the main assumption of our RG approach is the possibility of taking the continuum 
limit, which neglects correlation effects due to a large on-site attraction U . These are important for the Hubbard 
model at or near half-filling. Otherwise, our approach provides the evidence in a concrete model of the pmposed 
connection between charge-density-wave or phase separation instabilites and anisotropic superconductivity^!!. We 
have to bear in mind that, in our low-energy theory, the condensate wavefunction turns out to have the same sign 
near (tt, 0) and (0, it) but it must have necessarily less strength far from the saddle points, leading to a highly 
anisotropic gap. Specially sensible is our prediction of phase separation at low filling and weak coupling, when the 
Fermi sea is degenerating towards a pair of straight lines near t' — 0.5t. The interplaji-between phase separation and 
anisotropic superconductivity has been also considered in the t — J and U-V modelala as well as in the three-band 
Hubbard modelEl Phase separation near a Van Hove singularity has been discussed in a model with electron-phonon 
interaction in Ref. 

With regard to real systems, the realization of a charge-density- wave in the model |Biay bare a direct relation with 
the recent experimental observation of such instability in two-dimensional interfacescil. In those systems the atoms 
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at the surface are arranged in a triangular lattice, but the most important point is the present recognition that there 
is no significant nesting of the Fermi line accounting for the instability. The form of the dispersion relation has been 
evaluated in several instances, and it shows that the Fermi line is close to the saddle points at the boundary of the 
Brillouin ZoneE3. 

From the theoretical point of view, one of the most remarkable features of the phase diagram in Fig. [|is the existence 
of a triple point where the three different phases coalesce. The phase diagram is actually a map of the different ground 
states of the system. If we look at phase separation as an instability in which no symmetry breakdown takes place in 
a microscopic scale, opposite to what happens in the other two phases, we may conjecture that the triple point must 
correspond to a quantum critical point of the system. In fact, it has to be possible to crossover smoothly from the 
superconducting ground state to the charge-density-wave state through the phase separation region. The crossing of 
the two ground states is the distinctive feature of the quantum phase transtition. 

It is very suggestive the similarity that the phase diagram of Fig. Qjpears with that of the repulsive t — t 1 Hubbard 
model at the Van Hove filling obtained with analogous RG techniquesQ. In that case there is a large boundary in the 
phase space between antiferromagnetic and d-wave superconducting phases, up to a point where the leading instability 
turns out to be ferromagnetism. The existence of these three phases has been confirmed by numerical methods. Before 



clear signature of 
as well as mean- 



the RG analysis had been undertaken, quantum Monte Carlo computations had already provide' 
ferromagnetism at t' — 0.47£a. On the other-Jaand, recent quantum Monte Carlo computation: 
field computations with effective interactions! 2 ^, are supporting the existence of a d-wave superconducting phase at 
intermediate values of t'. These evidences prove the predictability of the RG approach to the present model, and they 
reassure us that the phases discussed above should be susceptible of being obtained by alternative numerical methods. 
The question of the mentioned quantum critical point is important enough for being tested by other computational 
techniques. A number of points, like the nature of the strong coupling phases or the thermodynamic properties of the 
model, should be also investigated to that effect. 

It is a pleasure to thank fruitful discussions with F. Guinea, S. Sorella and M. A. H. Vozmediano. This work has 
been partially supported by the CICYT Grant PB96-0875. 
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FIG. 1. Energy contour lines about the Fermi level, with the Fermi line passing by the saddle points A and B. 

FIG. 2. Different interaction terms arising from the flavor indexes A and B. 

FIG. 3. Second order diagrams renormalizing the interactions in the model, with electron lines carrying flavor index A or B 
appropriate to each case. 

FIG. 4. Phase diagram in the (t',U) plane showing the regions of charge-density- wave (CDW), s-wave superconductivity 
(SC-s) and phase separation (PS) instability. 
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